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Extensive first principles calculations are carried out to investigate gold-promoted TiO2(110) sur- 
faces in terms of structure optimizations, electronic structure analyses, ah initio thermodynamics 
calculations of surface phase diagrams, and ah initio molecular dynamics simulations. All com- 
putations rely on density functional theory in the generalized gradient approximation (PBE) and 
account for on-site Coulomb interactions via inclusion of a Hubbard correction, PBE+U, where U 
is computed from linear response theory. This approach is validated by investigating the interac- 
tion between TiO2(110) surfaces and typical probe species (H, H2O, CO). Relaxed structures and 
binding energies are compared to both data from the literature and plain PBE results, thus allowing 
the performance of the PBE+U approach for the specific purpose to be verified. The main focus 
of the study is on the properties of gold-promoted titania surfaces and their interactions with CO. 
Both PBE+U and PBE optimized structures of Au adatoms adsorbed on stoichiometric and re- 
duced Ti02 surfaces are computed, along with their electronic structure. The charge rearrangement 
induced by the adsorbates at the metal/oxide contact are also analyzed in detail and discussed. By 
performing PBE+U ah initio molecular dynamics simulations, it is demonstrated that the diffusion 
of Au adatoms on the stoichiometric surface is highly anisotropic. The metal atoms migrate either 
along the top of the bridging oxygen rows, or around the area between these rows, from one bridging 
position to the next along the [001] direction. No translational motion perpendicular to this direc- 
tion is observed. Approximate ah initio thermodynamics predicts that under 0-rich conditions, 
structures obtained by substituting a Tisc atom with an Au atom are thermodynamically stable 
over a wide range of temperatures and pressures that are relevant to applications in the realm of 
catalysis. Finally, it is shown that TiO2(110) surfaces containing positively charged Au ions activate 
molecular CO, whereas a single negatively charged Au"*^ species bound to an O vacancy only weakly 
interacts with CO. Despite this, the calculations predict that the reactivity of gold nanoparticles 
nucleated at O vacancies can be recovered for cluster sizes as small as Au2. 

PACS numbers: 68.43.Fg, 73.20.Hb, 68.47.Gh, 82.65.+r 



INTRODUCTION 



Titania, Ti02, is a metal oxide of both fundamental in- 
terest and technological importance^ .^. It is used in sev- 
eral key technologies including pigments, coatings, elec- 
tronic devices, implants, gas sensors, photochemical reac- 
tions, and catalysis ^^. One of the most important prop- 
erties of titania is that it can be easily reduced (see e.g. 
Ref. M for a concise presentation), strongly affecting its 
chemical properties in general and its reactivity in partic- 
ular—. One way to reduce the Ti02 surface is to remove 
surface oxygen atoms, thereby creating O vacancies. The 
removal of an O atom gives rise to two excess electrons 
and the appearance of new electronic states within the 
band gap at about 0.7 to 0.9 eV below the conduction 
band edge, thus creating an F-center— i^iiiii^. By this 
process, two substrate Ti^+ ions change formally to a 
Ti^+ oxidation state; see Ref. [l3l for recent literature and 
a detailed picture of the (de-) localization dynamics of 
the excess electrons. Alternatively, the Ti02 surface can 
be reduced by hydroxylation of surface O atoms via ad- 
sorption of hydrogen™— . The interaction of Ti02 with 
water is an important process which has to be taken into 
account, since it occurs easily, even in well-controlled 
UHV experiments. The adsorption of water on Ti02 has 
been investigated extensively, both experimentally and 



theoreticallj^iiSii^ — ; in particular see Refs. [TqIITqI for 
the most recent reviews of this literature. 

Most relevant to catalysis is the interaction of Au and 
CO with stoichiometric or reduced Ti02 surfaces and, 
in particular, the interaction of titania-supported gold 
particles with CO molecules. A detailed understanding 
of the process of CO adsorption is required to best com- 
prehend its wide variety of application^!^, such as CO 
oxidation at low temperature^^, the water gas shift re- 
action, and CO hydrogenation^ii^. In a recent paper—, 
we investigated the interaction of CO with the stoichio- 
metric TiO2(110) surface using a combination of density 
functional theory (DFT) and post Hartree-Fock meth- 
ods. For a single CO molecule in the (4x2) surface 
unit cell of our slab, we found that the upright position 
above the fivefold coordinated Ti sites, Tisc, remains the 
preferential adsorption geometry, even without enforc- 
ing symmetry. On the reduced titania surface, results 
from temperature-programmed desorption (TPD) exper- 
iments^*^ suggested that, at low coverages, CO adsorp- 
tion occurs at non- adjacent Tisc sites. These findings 
were supported by various calculations^^^^. However, 
earlier studies implicated bridge-bonded oxygen vacan- 
cies as adsorption sites for CO^i^, a conclusion corrob- 
orated by some theoretical investigations^!^ as well. 

The seminal work of Haruta and coworkers^ has 



shown that the low-temperature oxidation of molecular 
CO can be efficiently catalyzed by highly dispersed Au 
nanoparticles supported on Ti02 surfaces^—. It is now 
recognized that gold nanoclusters, prepared in different 
ways and supported on various metal oxides, are able to 
catalyze a number of reaction^*^, and that the size of 
the gold particles substantially affects the catalytic activ- 
ity. The gold clusters should be smaller than about 5 nm 
for high catalytic activity to occur, suggesting the key 
importance of metal/support interfacial interactions on 
a nanometer scale. Extensive studies of the Au/Ti02 sys- 
tem link the peculiar catalytic activity of gold nanopar- 
ticles on titania to several factors: high concentration of 
low-coordination site^*^, quantum size effects of two- 
layer Au islands^, active perimeter sites of the nanopar- 
ticle^, and charge transfer between the gold particles 
and the supporting oxide^*^. 

Over the past decade, DFT-based calculations have 
been extensively employed to study the interaction be- 
tween gold and the TiO2(110) surfaceii^Ii. Most of the 
existing theoretical studies provide information on stable 
adsorption sites of Au on the stoichiometric and reduced 
titania studies, while less effort has been devoted to the 
study of the 0-rich Au/TiO2(110) system^^^'^^ and the 
diffusion of Au adatoms on the stoichiometric and re- 
duced TiO2(110) surfaceii^i^. A wide variation in the 
lowest-energy positions of Au on titania are reported in 
the literature^, which can be explained in part by con- 
sidering that Au can diffuse rather easily on the stoichio- 
metric surface^. The potential energy surface (PES) 
of a single Au adatom deposited on the stoichiometric 
TiO2(110) surface or adsorbed into a surface O vacancy 
has been explored using static calculations^^. It has been 
shown that Au migration on the stoichiometric surface is 
two-dimensional, with a relatively flat profile. In the 
scenario where an Au atom is substituted for a surface 
Tisc site, it has been demonstrated that the Au atom is 
capable of weakening bonds of surface oxygens with the 
oxide^. 

Most of the density functional theory studies available 
in the literature dealing with defects and/or molecules 
adsorbed on titania substrates using reasonably sized su- 
percells make use of local (LDA) or semilocal (GGA) 
functionals. Despite widespread use, such functionals are 
known to often (but not always^) fail to predict quali- 
tatively correct electronic structures for reduced transi- 
tion metal oxides, due to the self-interaction error in- 
herent in the functionals. To partially correct for the 
self-interaction error, different computational methods 
can be used: perturbative many-body theories such 
as "GW — — , LDA plus dynamical mean field theory 
(DMFT)^, pseudo self-interaction-correction schemes 
(pSIC)2^, LDA plus Ull, and other methods that rely 
on hybrid functionals^^. Recently, GGA+U approaches 
have been applied with promising results in studies of 
intrinsic electron transport in Ti02 bulk^*^ and in the 
investigation of the charge (de-) localization dynamics in- 
duced by surface oxygen vacancies on the (110) surface 



of Ti02 in the rutile structure!^. Within the GGA+U 
approach, the electronic structure is partially corrected 
for the self-interaction error by adding a Hubbard term 
acting on the Ti-3^ orbitals. This approach has the ad- 
vantage of not adding much computational overhead to a 
standard GGA calculation in the plane wave / pseudopo- 
tential framework, thus enabling one to use fairly large 
supercells in order to allow structural relaxation to occur 
or to carry out ab initio dynamics. 

In this article, the PBE+U formalism is employed to 
investigate the structural and electronic properties of 
gold-supported TiO2(110) surface catalysts and their re- 
activity towards CO. As a first step, the PBE+U ap- 
proach is applied to study the interaction between the 
stoichiometric or reduced TiO2(110) surface and small 
probe species: H, II2O, and CO. The PBE+U structures 
and binding energies are compared to both the corre- 
sponding plain PEE results and reference data in the lit- 
erature in order to assess the applicability of the PBE+U 
approach for the present purpose. When it comes to gold 
on titania, our PBE+U ab initio molecular dynamics sim- 
ulations extend existing static relaxations and nudged 
elastic band mappings of the PES and demonstrate that 
Au adatoms diffuse in a highly directional manner on 
the stoichiometric surface. The metal atoms migrate 
easily, either along on top of the bridging oxygen rows 
or around the area between these rows, from one bridg- 
ing position to the next, along the [001] direction. The 
relative thermodynamic stability of different TiO2(110) 
structures is furthermore studied by employing the for- 
malism of approximate ab initio thermodynamics. Our 
calculations greatly extend existing studies and show that 
under 0-rich conditions, the thermodynamically most 
stable structure is the defective Au^Vtisc surface struc- 
ture (obtained by substituting a surface Tisc atom with 
an Au adatom), while under Ti-rich conditions, the Au 
adatoms are preferentially adsorbed at O vacancies. 

The remainder of the paper is organized as follows: In 
Section [III the model system, the methods and the com- 
putational details are summarized. The PBE+U formal- 
ism is validated in Section Hill bv studying the interaction 
of the TiO2(110) surfaces with a set of well-studied ad- 
sorbates such as H, H2O and CO. The main part of 
the paper is Section IIVI which presents novel insights 
into the structures, electronic properties, thermodynam- 
ics and dynamics of gold-promoted rutile TiO2(110) sur- 
faces and their interaction with CO. The concluding Sec- 
tion IVII summarizes the main results and puts them in a 
broader perspective. 



II. METHODS AND COMPUTATIONAL 
DETAILS 

The TiO2(110) surfaces were modeled by four 
O — Ti202 — O trilayer (4x2) supercell slabs separated 
by more than 10 A of vacuous space normal to the sur- 
face. The bottom of the slab was passivated with pseu- 



dohydrogen atoms of nuclear charge +4/3 and +2/3 in 
order to achieve well-converged results. This is our so- 
called "standard setup" which has been previously care- 
fully constructed^ by performing extensive tests on the 
convergence of surface energies as well as hydrogen and 
water adsorption energies, with respect to both the num- 
ber of relaxed outermost trilayers and the thickness of the 
slab itself (see tables and graphs in^ for detailed compar- 
isons). The system size employed in our calculations, 
corresponding to 208 atoms for the stoichiometric slab, 
belongs to the largest systems used so far in order to 
model the surface, in particular when it comes to per- 
forming ab initio molecular dynamics. In order to fur- 
ther check the convergence we optimized two five trilayer 
slabs (one with an empty surface O vacancy and one with 
an Au adatom at this vacancy) and confirmed that the 
resulting spin density and excess charge localization is 
the same as reported herein for our "standard setup" . 

The gradient-corrected Perdew-Burke-Ernzerhof func- 
tional (PBE)^^ was employed to describe semilocally the 
exchange-correlation effects. The spin-polarized Kohn- 
Sham equations were solved in the plane wave / pseu- 
dopotential framework using Vanderbilt's ultrasoft pseu- 
dopotentials^^ with a cutoff of 25 Ry using the F-point. 
The Ti pseudopotential was constructed from an ionic 
3d}4:S^ configuration and the 3s and 3p semicore elec- 
trons were treated as full valence states. It is well es- 
tablished that adding a Hubbard U term acting on the 
Ti-3d orbitals greatly improves the quality of LDA or 
GGAs in describing the electronic structure of both ox- 
idized and reduced titania surface o-^^i'^^i^^i^^' — . Follow- 
ing our previous work^^, we used a self-consistent linear 
response formalis m^^i^^ to compute the Hubbard term, 
which turns out to be U = 4.2 eV for this particular 
setup; the occupations of the d orbitals were calculated 
using atomic-like wave function projectors. It will not 
have escaped attention that our value of the U param- 
eter is larger than that recently derived by Mattioli et 
al^ (i.e. U = 3.25 eV). This can be attributed to the 
different (i-orbitals used as projectors for the integration 
of the (i-orbital occupation numbers. The U value of 
3.25 eV obtained in^^ was derived using the (i-orbital of 
the neutral Ti atom as a projector, whereas our value, 
U = 4.20 eV, is computed by using Ti+^ as a reference, 
which in our opinion more closely resembles the charge 
state of Ti in the TiO2(110) surface. At this point it 
should be noted that we also obtained U = 3.20 eV when 
using the (i-orbital of the neutral Ti atom as the projector 
instead, which is in agreement with the U value reported 
in^. Similarly, performing the calculations of FeO Pick- 
ett et al.^^ showed that the calculated Hubbard energy 
strongly depends on the choice of (i-orbital. Using the d- 
orbital of the neutral Fe atom they obtained U = 4.6 eV 
whereas using the Fe+^ dication yields a substantially 
larger value of U = 7.8 eV. This is a well-known and 
still poorly understood shortcoming of the U parameter 
derivation procedures that use atomic-like (i-orbitals as 
projectors for the integration of the occupations of the 



(i-orbitals in solids, which, in turn, is an input for the 
computation of the Hubbard correction to LDA/GGA 
density functionals^S. The static optimizations for the 
different TiO2(110) surface structures were carried out 
using the Quantum Espressc^^ code. All structures were 
relaxed by minimizing the atomic forces, where conver- 
gence was assumed to have been achieved when the max- 
imum component of the residual forces on the ions was 
less than 0.02 eV/A. Here, only the lowest trilayer atoms 
were constrained to their equilibrium positions while all 
other atoms were free to move during optimization. All 
ab initio molecular dynamics (AIMD) simulations^^ were 
carried out using the same spin-polarized PBE+U ap- 
proach, together with the Car-Parrinello propagation 
scheme, ^^ using a fictitious electron mass of 700 a.u. and 
a time step of 0.145 fs. Our in- house modified version of 
the CPMD^^ code was used for this purpose. 

The adsorption energy E^^^ per H atom on the stoi- 
chiometric TiO2(110) surface is computed from 
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where E^^^ ^ ^(^h) is the total energy of the slab satu- 
rated with A^H H adatoms; here E^^ ~ ^ ^ is the energy 
of the stoichiometric slab, which we take as a reference, 
and E^^l is the energy of a H2 molecule. When deal- 
ing with CO, H2O, and Au, the adsorption energies on 
stoichiometric and reduced TiO2(110) surfaces were cal- 
culated according to 
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where E^^l"^^, ^tot^^ ^^<^ ^t^t ^^^ ^^e total energies of 
the combined system, the (Au/)TiO2(110) surface in a 
certain oxidation state, and the isolated X adsorbate, re- 
spectively. The adsorption energies were calculated with 
and without inclusion of the Hubbard U term correction 
to the standard density functional, i.e. using the plain 
PEE and the PBE+U approaches. The 0-vacancy for- 
mation energy was calculated using 
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where E^^~^ and ^£^tot represent the total energy of 
the defective system and of the O atom, respectively. 
Because (semi) local functionals are known to overbind 
molecular O2, the total energy of the O atom was ad- 
justed in the manner of our previous work^. 

In order to analyze the thermodynamic stability of our 
different structures in the presence of H adatoms, we em- 
ploy the formalism of approximate ab initio thermody- 
namics^ ^i^ by assuming that the surfaces can exchange 
H atoms with a surrounding gas phase. Assuming ther- 
modynamic equilibrium, the most stable surface compo- 
sition at a given temperature T and pressure p is given 
by the minimum of the surface Gibbs free energy. Since 



we are only interested in the relative stabilities of sur- 
face structures, we directly compute the differences in 
the surface Gibbs free energies AGads{T^p) between the 
defective and the ideal surface according to 
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where A is the surface area, ATVh is the difference in 
the number of H atoms between the two surfaces, and 
Mh(^,p) is the chemical potential representing the Gibbs 
free energy of the gas phase with which the H atoms are 
exchanged. Assuming that all differences in entropy and 
volume contributions in AG^dsiT^p) are negligible, the 
Gibbs free energies are approximated by their respective 
total energies of our DFT slab calculations as usua^ ^iSl. 
The upper bound for the chemical potential jiYiiT^p) is 
given by the total energy of its most stable elemental 
phaseS^, that is, molecular hydrogen {^Ef^^). This upper 
bound is taken as the zero of our energy scale by using 

In a similar way, the effect of temperature and pres- 
sure on the relative stability of the Au/TiO2(110) sur- 
face structures is studied by employing the formalism of 
approximate ah initio thermodynamics^ -i^. The free 
energy of formation of the Au/TiO2(110) surface struc- 
tures AadsG^(^7P) is assumed to depend on the tempera- 
ture and pressure only via the oxygen chemical potential 
l^o{T,p) given by 



/io(T,p)=/io(T,/) + -feTln 



(5) 



Equation (|5]) represents the thermodynamics reservoir 
of the O2 environment that is in contact with the sur- 
face under consideration. The free energy differences will 
be calculated as a function of AjioiT^P) = fio(T^p) — 
/xo(^ = K,p^), corresponding to changes of the oxygen 
chemical potential with respect to a zero reference state. 
The latter is set to the total energy of the O atom at 
T = K, /io(T = K,pO) = 1/2E2^^ = 0. Assuming 
thermodynamic equilibrium of the surfaces with an O2 
gas phase, the chemical potential can be converted into 
a pressure scale for different temperatures by using ex- 
perimental thermochemical reference data or by applying 
the ideal gas equation22rJ£i. Vibrational and rotational 
entropic contributions to (iq (T, p) are include d by means 
of thermodynamic tables as described in Ref. [99|. Under 
these assumptions and neglecting entropic contributions 
of the solids involved, the free energy of formation as a 
function of pressure and temperature assumes the expres- 
sion 
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where E^q- is the energy of a formula unit of the Ti02 
bulk phase. The quantities Nq and A/'^J- represent the 
number of O or Ti vacancies that are present in the struc- 
ture under consideration. Therefore, the energy cost for 
the formation of surface defects is taken into account in 
Eq. (j6]) via the chemical potential of O atoms and of bulk 
Ti02. Finally, the chemical potential of Au, /iAu, is set to 
be the total energy per atom of the bulk Au crystal. The 
upper bound for the chemical potential jiq is given by the 
total energy of its most stable elemental phase, that is, 
molecular oxygen (^E^^^). This upper bound is used to 
define AjioiT, P) = jjlq — \E^^^. A lower bound for A/io 
is given by minus half of the formation energy of bulk 



Ti02, i.e. E/'"^ = ^buik + ^to? - ^bulk^ for which we 
have taken the theoretical value of 4.8 eV from our PBE 
calculations; here E^^^^ and ^^^uik ^^^ ^^^ energies of one 
bulk unit cell of Ti02 and metallic Ti, respectivelySSiiffi. 
Finally, the bonding charge density has been evaluated 
using the expression 



Apif) = psub+X - (psub + Px) 



(7) 



A^Ti[^Ti02-2M0(T,p)]-/iAu] 



where the p's are the respective valence electronic charge 
densities at position r in space. 



III. SMALL MOLECULES ON TITANIA: 
REFERENCE CALCULATIONS USING PBE-hU 

Computationally efficient implementations of DFT 
based on local/semilocal LDA/GGA density functionals 
predict rather delocalized defect levels for excess elec- 
trons in the case of reduced transition metal oxides in 
general, and for titania in particular. Thus, more sophis- 
ticated techniques such as hybrid functionals or GGA+U 
approaches are necessary to properly account for the 
strong correlation effects of these (i-electrons, resulting in 
localization of the excess charge on 3<i-orbitals of reduced 
Ti atoms. However, before using PBE+U to investigate 
the properties of gold on titania, it is necessary to confirm 
that this approach does not destroy the agreement be- 
tween previously reported plain PBE results (mainly op- 
timized structures and relative energies) and experimen- 
tal observations. The purpose of this section is therefore 
twofold: first, to validate the PBE+U approach using a 
test set that probes the physics and chemistry of titania 
surfaces interacting with adspecies relevant to heteroge- 
neous catalysis; and, second, to check how PBE+U per- 
forms compared to plain PBE investigations for the very 
same systems. Thus, the interaction of H, H2O, and CO 
with the TiO2(110) substrate is investigated in the fol- 
lowing section using the PBE+U approach where a large 
set of both experimental and previous theoretical data 
are available. 

Many theoretical and experimental studies have been 
devoted to understanding the interaction between hydro- 
gen atomsSii^i^ii22 or water^^J^J^^ and the TiO2(110) 
surface. In experiments, hydrogen atoms, when adsorbed 



on TiO2(110), stick to the bridging oxygens and a maxi- 
mum surface saturation limit of ~ 0.7 ML is observed^-. 
These findings have been confirmed on purely theoretical 
grounds within the framework of standard GGA calcula- 
tions using the PBE functional^. In-depth ab initio ther- 
modynamics considerations reveal a maximum saturation 
level of hydrogen on this oxide surface of about 60-70%, 
in excellent agreement with the above-mentioned experi- 
mental observations. Adsorption of hydrogen on the sto- 
ichiometric surface results in its reduction by introducing 
one electron per adsorbed H atom into the substrate. 

We first consider the interaction of H with the stoichio- 
metric TiO2(110) surface. The hydroxylated surface has 
been investigated for a wide range of H coverages using 
both PBE-FU and PBE. In agreement with both plain 
PBE calculations and experimental results, we find that 
H atoms preferentially adsorb on top of surface Ob atoms 
of the stoichiometric titania surface, leading to the for- 
mation of OH groups with 0-H bond lengths of ~ 1 A. In 
addition, our PBE+U calculations predict that one elec- 
tron per H atom adsorbed on the surface is transferred 
to the substrate, leading to the formation of Ti^+ ions 
(see Table HI). The newly- formed OH groups are found 
to be tilted by about 20-50° in opposite [110] directions, 
as previously reported. The PBE and PBE+U values 
of the adsorption energies E^^^ per H atom at different 
coverages are reported in Table HI The ground state con- 
figuration of the fully hydroxylated titania surface yields 
(2x1) symmetry, but at room temperature OH groups 
will be fully disordered with respect to their axes be- 
cause of the tiny barrier that must be overcome in order 
to flip their orientation. 

At all coverages, the PBE+U adsorption energies are 
found to be significantly lower, by about —0.3 to —0.6 eV, 
compared to the PBE data in Table HI However, the same 
stability trend obtained by employing the PBE+U func- 
tional is obtained using the standard PBE functional, 
i.e. a decrease of the adsorption energy with increasing 
coverage. Nevertheless, even at full monolayer coverage, 
the adsorption energy per atom is still significant and 
negative when using PBE+U, whereas it is close to zero 
according to PBE. Clearly, such total energy considera- 
tions need to be supplemented with ab initio thermody- 
namics in order to check the thermodynamic stability of 
the surface at different saturation levels. The results are 
summarized in terms of the surface free energy diagram 
in Fig. [H Both PBE+U and PBE predict that the fully 
hydrogenated surface will be thermodynamically unsta- 
ble at all accessible hydrogen partial pressures. PBE+U 
and PBE also agree in that saturation of this surface 
by hydrogen is reached at a coverage on the order of 
70%. Thus, the previous PBE predictions^ and agree- 
ment with experimental observations^^ upon hydroxy- 
lating TiO2(110) are qualitatively confirmed, (for rea- 
soning, see the detailed conceptual discussion in Ref. |9|)- 
However, significant strengthening of the adsorption and 
shifts in the surface phase diagram are observed when 
using PBE+U, thus localizing the excess charges upon 



TABLE I: Adsorption energies E^ds (in eV) per H atom for 
adsorption of hydrogen at different coverages as indicated. In 
the second column the number of reduced Ti"^^ ions which 
are present in the substrate is reported. 
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introducing the Hubbard correction. 
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FIG. 1: Free energy AGads{T,p) for H adsorption on 
TiO2(110) stoichiometric surfaces with different hydrogen 
coverages as a function of the hydrogen chemical potential 
A/iH- Conversion to hydrogen partial pressures pu2 (upper 
axis) has been carried out at T = 400 K (see text). The 
red and blue lines represent the PBE+U and PBE results, 
respectively. 



It is well-known that the investigated surface contains 
a significant number of oxygen vacancies {^ 5 %), not 
just in ill-defined industrial situations, but even under 
well-controlled experimental conditions^ Therefore, the 
investigation of the adsorption of water molecules on the 
reduced TiO2(110) surface is of great importance in the 
frame of addressing the fundamental issue of dissocia- 
tive versus molecular adsorption mode^ i-^^i-^^i^^ . We thus 
consider the interaction between H2O and the Ti02 sur- 
face. We have performed systematic PBE+U and PBE 
calculations of H2O adsorbed on the reduced Ti02 sur- 
face, considering a titania surface containing a single Vq 
vacancy; this greatly extends our recent comprehensive 
work^ concerning water on the stoichiometric surface us- 
ing plain PBE. 

The specific adsorption configurations for a water 
molecule on the reduced surface, which we consider in 



this investigation, are compiled in Fig. [21 and the corre- 
sponding adsorption energies are cohected in Table |TT1 
The Tisc surface sites are labeled in relation to the 
Vo vacancy, where site TiO denotes a nearest-neighbor 
Tisc atom and sites Til and Ti2 are the second and third 
nearest-neighbour Tisc sites parallel to the Ob row con- 
taining the Vo site; see panel (a) of Fig. [2j In agreement 
with experimental data, we find that on the reduced sur- 
face, H2O molecules prefer to dissociatively adsorb at Vo 
vacancy sites, leading to a configuration with two surface 
OH groups as shown in panel (b) of Fig. O We therefore 
end up with a stoichiometric titania surface with two H 
atoms adsorbed on two surface Ob atoms. Once the water 
molecule dissociates at the oxygen defect through proton 
transfer to an adjacent Ob atom, the PBE+U (PBE) 
adsorption energy is —1.61 eV (—1.18 eV). A projected 
PDOS analysis reveals that upon dissociation of water at 
the Vo site, two second-layer Ti^+ ions are present in the 
substrate. This value of Eads is ~ 0.4 eV lower than the 
corresponding plain PBE value and previously reported 
valuesr^i^ii^ which are between —0.94 and —1.1 eV. 
The energy value deduced from a water desorption peak 
at 520 K in TDS experimentsii^ using the simple Red- 
head formula^ is about —1.4 eV, which is between the 
PBE+U and PBE values. However, the estimation of 
desorption energies using the Redhead formula can be 
biased by as much as 25%, which implies that both val- 
ues must be considered to be consistent with experiment. 
We now turn our attention to the adsorption and disso- 
ciation of water at Tisc sites next to Vo • We anticipate 
that, as observed in the case of H2O dissociatively ad- 
sorbed at the Vo site, the interaction between water and 
the Tisc sites does not further reduce the metal oxide 
support. All PBE+U calculations predict the presence 
of two reduced Ti^+ ions before and after the adsorption 
of H2O at Tisc sites. Water molecules can be adsorbed ei- 
ther dissociatively (labeled as "D" ) or molecularly ( "M" ) 
at the various Tisc sites (i.e. TiO, Til, and Ti2) next to 
Vo (see Fig.[2j). When H2O dissociates at a Tisc site, the 
resulting configuration contains an OH group bonded to 
a Tisc atom and an H atom bonded to a nearest-neighbor 
O atom of the Tisc in the [110] direction. We have con- 
sidered two different topologies: first, where the H atom 
coming from the dissociated water molecule binds to an 
O atom belonging to the Ob row in which the Vo vacancy 
site is present (labeled configuration "A"); and second, 
with the H atom bonded to O atoms belonging to an 
adjacent Ob row parallel to the Ob row that hosts Vo 
(configuration "B"). See Table ITTl for the corresponding 
adsorption energies. 

Our PBE+U and PBE calculations suggest that, in the 
presence of an oxygen vacancy, water molecules adsorb 
dissociatively or molecularly at Tisc sites with £^ads in 
the range of —0.45 to —0.92 eV. Once the water molecule 
dissociates at a Tisc site, it forms a pair of terminal hy- 
droxy Is. Therefore, in the dissociative case, we always 
end up with an OH group bonded to a Tisc site and a 
protonic H atom transferred to an adjacent Ob atom. 



TABLE U: Adsorption energies ^ads (in eV) in the case of 
adsorption of water on the reduced surface for PBE+U and 
plain PBE reported in parentheses; see Fig. [2] for labeling. 



configuration 


site TiO 


site Til 


site Ti2 


DA 


2.27 ( 1.58) 


-0.71 ( -0.81) 


-0.45 ( -0.48) 


MA 


-0.77 ( -0.79) 


-0.73 ( -0.76) 


-0.72 ( -0.72) 


DB 


-0.87 ( -0.89) 


-0.92 ( -0.95) 


-0.89 ( -0.90) 


MB 


-0.77 ( -0.80) 


-0.83 ( -0.81) 


-0.84 ( -0.82) 


H20@Vo 


-1.61 ( -1.18) 







the Ob atom belonging either to the Ob row that con- 
tains Vo, or to the adjacent one that contains no oxygen 
vacancy. As shown in Table [III the most stable dissocia- 
tive configurations are those with the protonic H atom 
transferred to an Ob atom of an Ob row parallel to the 
row featuring the Vo site (configurations denoted "DB" 
in Table inland Fig. [2]). The PBE+U/PBE calculations 
show that protons prefer to bind to Ob atoms belonging 
to an Ob row in absence of Vo sites, with £^ads ranging 
from -0.87 to -0.95 eV. 

Otherwise, if we consider the case where, after water 
dissociation at a Tisc site, a proton transfers to an Ob 
atom belonging to an Ob row that includes a Vo site 
(configurations denoted "DA" in Table HIl and Fig. [2j), 
our calculations provide adsorption energies in the range 
of —0.45 to —0.71 eV; the same trend is observed with 
or without the inclusion of a Hubbard U term in the 
calculations. We note that the structure with an OH 
group adsorbed at a Tisc site and an H placed right at 
the Vo vacancy is unstable; the adsorption energy of this 
configuration £^ads is positive by 2.27 eV. 

Our results show that on the reduced TiO2(110) sur- 
face, water prefers to adsorb dissociatively onto Tisc sites. 
Once water molecules are molecularly adsorbed on the re- 
duced surface, the binding energy of the surface is about 
0.1 eV higher than it is when water is dissociated on the 
same substrate. These findings are in agreement with 
previous studies of H2O interaction with the stoichio- 
metric TiO2(110) surface. However, as shown based on 
carefully converged calculations^, molecular and dissoci- 
ated configurations become essentially energetically de- 
generate at very low coverages, which explains why some 
studies favor molecular adsorption whereas others yield 
the dissociated state as the lowest energy configuration in 
this regime^^. Most experimental works, except a recent 
Qj^gioe claiming mixed adsorption, indicate molecular ad- 
sorption only. This gives rise to the well-known discrep- 
ancy between theoretical predictions and experimental 
results in providing a consistent and comprehensive pic- 
ture of water adsorption on titanium dioxide surfaces, in 
particular at low coverages. 

Last but not least, we focus on the interaction between 
CO and the TiO2(110) surface. It is well known that, on 
the stoichiometric TiO2(110) surface, the CO molecule 
adsorbs onto Tisc sites, thus forming Ti-C bonds. Be- 
cause van der Waals dispersion interactions and non- 




at site 2 



FIG. 2: Ball and stick models of relevant configurations 
(see text) for an H2O molecule adsorbed either molecular ly 
("M") or dissociatively ("D") in two configurations ("A" and 
"B") on TiO2(110) surfaces (top view) at Vo, TiO, Til, and 
Ti2 sites (see panel (a) for site labeling) obtained using the 
PBE+U approach. Red, blue, violet, and yellow spheres are 
substrate O, Ti, water O, and H atoms, respectively, and the 



TABLE III: Adsorption energies E^^ds (in eV) of a CO 
molecule adsorbed on the reduced rutile TiO2(110) surface 
at different sites, labeled according to Fig. [S] In the second 
column, the number of reduced Ti^^ ions present in the sub- 
strate is reported. 

Configuration ^ads (PBE+U) Ti^+ ^ads(PBE) 



CO@Vo 


-0.32 


2 


-0.29 


CO@TiO 


-0.22 


2 


-0.28 


CO@Til 


-0.29 


2 


-0.30 


C0@Ti2 


-0.33 


2 


-0.29 



local electron correlations significantly influence this type 
of bonding, we have previously investigated the interac- 
tion of CO with the stoichiometric TiO2(110) surface us- 
ing a combination of DFT and post Hartree-Fock ("SCS- 
MP2") methods^^. The CO binding energy has been 
found to vary significantly with coverage and increases 
upon reaching the saturation limit. The SCS-MP2 Eads 
computed values are —0.20 eV for the full saturated sur- 
face and —0.36 eV for a single CO molecule adsorbed 
on the surface. The PEE adsorption energy of a single 
CO molecule, -0.32 eV, is close to the SCS-MP2 value of 
—0.36 eV, and both energies are in accord with the exper- 
imental value obtained by means of thermal desorption 
spectroscopy. This demonstrates that PEE as such is 
able to describe the interaction of CO with the ideal (110) 
rutile surface in the limit of low coverages, which remains 
unaltered when using PEE+U, which gives —0.31 eV for 
the adsorption energy. 




FIG. 3: Ball and stick model of the TiO2(110) surface (top 
view), with site labeling. Red and blue spheres are O and Ti 
atoms, respectively, and the oxygen vacancy site, Vo, within 
the bridging oxygen row, Ob, is highlighted using a green 
sphere. 



We now turn our attention to CO adsorbed onto the 
O vacancy site, Vq, as well as on the fivefold coordi- 
nated surface sites, Tisc, at various distances from the 
vacancy. Note that the Tisc site belongs to the Ti row 
next to the bridging Ob row containing Vo (see Fig. [3|). 
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FIG. 4: Electronic structure analyses (based on the PBE+U approach) of molecular CO adsorbed on the reduced TiO2(110) 
surface at the Ti2 site with respect to the oxygen vacancy, Vo, located in the row of bridging oxygens, Ob (see Fig. [3] and 
Table [nil). The left panel represents the bonding charge Sp(^z) integrated in planes perpendicular to the surface and plotted as 
a function of the height from the surface. The central panel displays the bonding charge Ap{r) at an isovalue of ±0.06 |e|/A^ 
where electron accumulation and depletion are represented by red and blue areas, respectively. The right panel shows the total 
DOS and atom-resolved projected DOS (PDOS) as indicated; here, energy values are with respect to the Fermi level, which is 
marked by a solid vertical line. 



In each case, the CO molecule is placed perpendicular 
to the substrate, with the carbon atom pointing toward 
the surface. All the structures are fully relaxed accord- 
ing to our aforementioned convergence criterion. The 
computed adsorption energies (Eads) are compiled in Ta- 
ble [Till Interestingly, we do not see a significant variation 
in the adsorption energies computed for different struc- 
tures. The PBE+U (PEE) adsorption energy for the 
CO molecule adsorbed at the surface Ob vacancy site 
is -0.32 eV (-0.29 eV). CO adsorption at sites Til and 
Ti2 results in adsorption energies of —0.29 eV (—0.30 eV) 
and —0.33 eV (—0.29 eV), respectively (see Fig.[3]for site 
labeling scheme). The computed values of Eads are in 
qualitative agreement with previous studiesiSiiifl^. We 
note that these values are close to those obtained when 
CO is adsorbed at the Vo vacancy site (—0.32 eV) and on 
the stoichiometric surface^, i.e. —0.31 eV (—0.32 eV) for 
PBE+U (PBE). When adsorbed at site TiO, the adsorp- 
tion energy of the CO molecule, computed with PBE+U, 
results in a distinctly higher value for Eads of —0.22 eV 
in agreement with Ref. 108, while the corresponding 
PBE value is —0.28 eV, comparable to the adsorption at 
sites Til and Ti2. Our PBE+U results therefore confirm 
previous findings that CO molecules weakly interact with 
the reduced TiO2(110) oxide surface, Eads being of the 
order of about -0.3 eV (see Refs. 39T07T08). The cal- 
culations indicate that CO adsorbs at both Vo vacancies 
and Tisc sites, but while PBE calculations give similar 
energy values for CO adsorption at Vo vacancies and at 
TiO, Til, and Ti2 sites, the inclusion of a Hubbard U term 
suggests that the adsorption of CO at TiO sites, the sites 
facing the Vo vacancy, is discouraged. In this case, the 
adsorption energy Eads is found to be ~ 0.1 eV higher 
when compared to the adsorption energy values at sites 
Til and Ti2. Upon CO adsorption on the reduced sur- 



face, the charge redistribution that results from attaching 
the molecule does not further reduce the oxide support: 
the PBE+U calculations yield two Ti^+ ions before and 
after adsorption, which results in an insignificant change 
in the adsorption energies when switching from PBE to 
PBE+U calculations. This behavior is confirmed by the 
computed electronic density of states (DOS). In Fig. 21 
we depict the electronic DOS and the bonding charge 
density Ap of the structure where a CO molecule is at- 
tached at the Tisc site (labeled site Ti2 in Fig. [3]). The 
electronic DOS features a peak in the band gap below the 
Fermi level. The projected DOS (PDOS) analysis reveals 
that this filled gap state is related to the charge localized 
on two Ti-3d orbitals that yield the two reduced Ti^+ 
sites. We conclude therefore that the binding of a CO 
molecule does not induce a significant charge rearrange- 
ment at the CO/oxide contact site. However, as we will 
show, the catalytic activity of the TiO2(110) substrate 
for efficient CO oxidation is improved by supported and 
dispersed Au adatoms on this substrate. 



In conclusion, this detailed assessment convincingly 
demonstrates that although similar trends are observed 
in the adsorption energies with or without the inclusion 
of a Hubbard U correction, the PBE+U method is seen 
to significantly improve the description of the electronic 
structure whenever reduction occurs. In particular, the 
localization of excess charge on the titania substrate in- 
duced by O vacancies (F-centers) or upon H atom ad- 
sorption on bridging O atoms (hydroxylation), is cor- 
rectly predicted by the PBE+U approach. Clearly, an 
adequate description of the electronic structure of such 
TiO2(110) surfaces is crucial when dealing with metal- 
promoted oxide surfaces in the realm of catalysis. 



IV. GOLD-PROMOTED TITANIA: 

ELECTRONIC STRUCTURE, DYNAMICS, AND 

THERMODYNAMICS FROM PBE+U 

A. Au adatom adsorption on the stoichiometric 
TiO2(110) surface 

Having shown that the PBE+U formahsm performs 
weh for a set of reference calculations on the adsorp- 
tion of H, H2O, and CO on stoichiometric and re- 
duced TiO2(110) surfaces, and that it significantly im- 
proves the description of reduced titania surfaces, we 
now progress to investigating the interactions between 
gold and TiO2(110) surfaces. The adsorption or substitu- 
tion of gold induces strong charge rearrangements at the 
Au/oxide contact, which affects the electronic structure. 
Of particular interest in the realms of metal/support in- 
teractions and heterogeneous catalysis is the oxidation 
state of Au adatoms, which is determined by the site 
where the metal atom is adsorbed, as well as by the 
stoichiometry of the supporting oxide. As in previous 
GGA studiesii^i^i^, two stable adsorption sites of a 
single Au adatom on the stoichiometric TiO2(110) sur- 
face have been identified, the two structures differing by 
only ~ 0.1 eV in energy. The most stable adsorption site 
for an Au adatom deposited on this titania surface is a 
bridge site between an Ob and a Tisc atom as depicted in 
the central panel of Fig. [5] (B). The computed PBE+U 
(PBE) adsorption energy and the Au-O / Au-Ti bond 
lengths are -0.58 eV (-0.41 eV) and 2.30 (2.39) / 2.79 A 
(2.88 A ) respectively, in agreement with previous stud- 
ie^^^i^i^ based on standard GGA calculations; see 
Table HV] for a summary. 

The bonding charge density analysis reveals that 
0.11 |e| are transferred from the metal atom to the oxide 
substrate, thus indicating a very weak oxidation of Au. 
The excess charge in the substrate is mostly localized 
around the Ob bonded to the Au adatom. This value of 
the charge transfer has been obtained by integrating the 
bonding charge density on planes parallel to the surface 
from the center of the vacuum region to the center of the 
O- Au-Ti bond (see left panel of Fig. [5] (B)). As demon- 
strated by the PDOS analysis shown in the right panel 
of Fig. [5] (B), in this configuration all the Ti ions belong- 
ing to the substrate preserve their formal oxidation state 
Ti4+. 

The second identified stable site is a top site (see cen- 
tral panel of Fig. [5] (C)), where the Au adatom is ad- 
sorbed on top of an Ob atom. The corresponding PBE+U 
(PBE) adsorption energy and the Au-O bond length are 
-0.48 eV (-0.29 eV) and 2.00 A (2.16 A), akin to previ- 
ous GGA studie^^. Again, a net charge transfer from 
metal to surface, leading to a positively charged Au'^^ 
ion, is observed. In this case, however, the magnitude 
of the charge transfer, 0.35 |e|, is more significant, i.e. 
three times larger than in the previous case. The excess 
charge in the substrate is now mostly localized around 
the surface Ob atom bound to the Au'^^ and a second- 



TABLE IV: Adsorption energies E^^ds (in eV) of Au and CO 
species on stoichiometric and reduced rutile TiO2(110) sur- 
faces and on the Au/Ti02 metal/support system; see text for 
labeling. In the second column, the number of reduced Ti^^ 
ions present in the substrate is reported. 





^ads (PBE+U) 


Xi3+ 


^ads(PBE) 


Au@0 (bridge) 

Au@0(top) 

Au@Vo 

Au2@Vo 

Au@VTi5c 


-0.58 
-0.48 
-1.54 
-1.17 
-6.38 



1 
1 
2 



-0.41 
-0.29 
-1.57 
-1.19 
-6.21 


CO@Au@0 (bridge) 

CO@Au@0(top) 

CO-Au@Vo 

C0-Au2@Vo 


-2.27 
-2.66 
-0.41 
-1.22 


1 
1 
1 
2 


-2.17 
-2.33 
-0.32 
-1.00 



layer Ti ion which reduces Ti'^+ -^ Ti^+ (see Fig.[5](C)). 
The reduced Ti^+ ion is located at a site adjacent to the 
Au'^^ adatom in the second subsurface layer under the 
Tisc. These findings are corroborated by the computed 
DOS plotted in Fig.[5](C), which displays two features in 
the band gap. The projected DOS analysis reveals that 
the filled state below the Fermi level and closest to the 
valence band results from the charge transferred from the 
metal to the substrate being localized on a second-layer 
reduced Ti^+ atom. The unoccupied level closest to the 
conduction band is instead related to the Au-O bonding. 
The 65 levels of Au are partially empty and are located 
above the Fermi level, leading to the Au oxidation. 



In summary, the PBE+U calculations predict two 
lowest-energy configurations for Au adsorption on the 
stoichiometric TiO2(110) surface: a bridge site with the 
Au adatom adsorbed between Ob and Tisc atoms and a 
top site with the Au adatom adsorbed on top of an Ob 
atom. Once adsorbed at the bridge site, a very weak 
oxidation of the Au adatom is observed. On the other 
hand, the adsorption process of Au on top of an Ob atom 
induces a net charge transfer from the adsorbate to the 
substrate, leading to the formation of a distinctly posi- 
tively charged Au'^^ species where about a third of an 
electron is transferred from the metal atom to the ox- 
ide substrate. A qualitatively similar scenario has been 
observed in recent studies of the related Au/Ce02 sys- 
ten>i22 liii. However, unlike titania, with ceria the charge 
transfer involved in the adsorption of Au on the stoichio- 
metric oxide surface always leads to the reduction of a 
substrate Ce ion. In addition, the fact that the excess 
charge (5—, stemming from Au in the present case, local- 
izes on a second-layer Ti ion is in line with our recent 
findings on reduced titania surfaced. 
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FIG. 5: Electronic structure analyses (based on the PBE+U approach) of an Au adatom (B) supported by the stoichiometric 
TiO2(110) surface in the bridge position, (C) supported by the stoichiometric TiO2(110) surface in the top position, (D) 
adsorbed on a surface Vo vacancy on the TiO2(110) surface, and (E) substituting a surface Tisc atom in the presence of a 
surface Vo vacancy in a bridging position on the TiO2(110) surface (see Fig.[9|). Left panels represent the bonding charge Sp(z) 
integrated over planes perpendicular to the surface and plotted as a function of the height from the surface. Central panels 
display the bonding charge Ap{f) at an isovalue of ±0.06 |e|/A^ where electron accumulation and depletion are represented 
by red and blue areas, respectively. Right panels show the total DOS and atom resolved projected DOS (PDOS) as indicated, 
where energy values are with respect to the Fermi level, which is marked by a solid vertical line. 
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B. Au adatom diffusion on the stoichiometric 
TiO2(110) surface 

Several previous theoretical studie^i^i^ explored the 
PES of a single Au adatom deposited on the stoichiomet- 
ric TiO2(110) surface or adsorbed onto a Vo vacancy site 
using static calculations, including nudged elastic band 
(NEB) mapping^i^. A key finding of these investiga- 
tions is that the PES for Au migration is quite fiat, with 
low energy barriers. This indicates that Au might dif- 
fuse rather easily on the stoichiometric surface. The esti- 
mated values of the energy barrier^ i^^i^^ agree with the 
experimental observation of facile Au diffusion on the ox- 
ide surface even at temperatures as low as 140 K, as well 
as the estimates for the binding energy of 0.5 eV and 
small migration barriers of 0.07 eV (see Refs. [12|15Q[ ). 

Inspired by these findings, we decided to perform ex- 
phcit dynamics using unconstrained ab initio molecular 
dynamics^ in order to reveal the mechanism of diffusion 
of an Au adatom on the stoichiometric TiO2(110) sur- 
face. In order to probe the dynamics more efficiently, the 
temperature of the simulations was set to T = 900 K us- 
ing the Car-Parrinello scheme^^ to propagate the system 
consistently using the PBE+U functional. The selected 
temperature is far above ambient yet sufficiently low so 
as to not decompose the surface. Thus, the phonon dy- 
namics is accelerated and the sampling of the PES is en- 
hanced on the picosecond AIMD time scale. As starting 
configurations for the AIMD simulations, we employed 
one structure in which the Au adatom is adsorbed on 
top of an Ob atom (see Fig. [5] (C) and Table [IVj) and a 
second structure where it is adsorbed in a bridge position 
between an Ob atom and a Tisc atom (see Fig.[5](B) and 
Table [IV|) . After equilibrating the structures at 300 K for 
several picoseconds, the system was heated to the target 
temperature of 900 K for the present analysis. 

Let us first consider the scenario with the Au adatom 
adsorbed on top of a surface Ob atom, labeled as site 01 
in Fig. [6l where the diffusion path of the Au adatom on 
the stoichiometric TiO2(110) surface is visualized. Dur- 
ing the simulation, the Au adatom diffuses in the [001] 
direction along the row of bridging oxygen atoms, Ob- 
Adatom diffusion is mediated by the Au atom hopping 
between nearest-neighbor oxygen atoms. As demon- 
strated in Fig. [6l the Au adatom is originally bonded to 
the surface bridging 01 atom with an Au-01 bond length 
of ~ 2 A (red line); it diffuses along the Ob row and after 
~ 0.75 ps reaches a configuration where it is equidistant 
between the 01 and 02 atoms. Then it jumps on top of 
the row's next atom, the 02 site, where the Au-02 bond 
length is ~ 2 A (green line). The Au diffusion proceeds 
along the Ob row and at about 2 ps the Au adatom is 
shared between 02 and the next site, 03 (blue line), until 
it jumps on top of 03 forming a bond of ~ 2 A. 

The charge localization and charge hopping dynamics 
along the adatom migration path is monitored by com- 
puting, as a function of time, the occupation matrix of 
each Ti d-a and (1-/3 spin orbital along the trajectory 



(same analyis as in^^). As shown in Fig. [8] the excess 
charge donated by the Au atom to the substrate is ini- 
tially localized on the second-layer Ti2 site but transfers 
from there to site Ti3 on the sub-picosecond time scale 
(at t ~ 0.75 ps). As seen by comparing Fig. [8] to Fig. [6l 
one observes that t ~ 0.75 ps corresponds exactly to the 
jump of Au from the 01 to the 02 site. A qualitatively 
similar scenario happens at about 2 ps which corresponds 
to the next hopping event of the Au atom from site 02 to 
03. Thus, the motion of the surface gold adatom along 
the row of bridging oxygen atoms, Ob, appears to be fully 
correlated with the localization and hopping dynamics of 
the excess charge injected into the oxide support in the 
second layer of Ti atoms. 

This dynamical scenario is distinctly different from 
what has been found recently for the excess charge in- 
duced by oxygen vacancies Vo in the bridging row on 
the same substrate—. In the presence of gold adatoms, 
the present simulations suggest a more localized configu- 
ration for the excess electron. This localized electron ap- 
pears to preferentially populate sites in the vicinity of the 
Au atom and to closely follow the motion of the oxidized 
adatom. Another interesting phenomenon observed dur- 
ing the simulation is the absence of excess charge on the 
substrate at about 1.5 ps. Computing the spin density 
close to t = 1.5 ps we observe that the charge localized at 
second-layer Ti sites disappears from the substrate and 
goes to the Au adatom, where it sits for a fraction of a pi- 
cosecond before returning to the substrate and occupying 
the Ti3 site. As observed in our previous work^^ the ex- 
cess charge populating specific second-layer Ti sites and 
coming from the Au adatom adsorbed on top of Ob atoms 
migrates easily by phonon-assisted (thermally activated) 
hopping to other Ti sites. 

Next we consider the situation where the Au adatom 
is initially adsorbed at a bridge position between an Ob 
atom and a Tisc atom (see Fig. [7j). There it forms two 
bonds with the 01 and Til atoms, with Au-01 and Au- 
Til bond lengths of - 2.3 A (red line) and - 2.8 A, 
respectively. Also in this case, the gold atom originally 
bonded to the 01 and Til atoms is found to diffuse ex- 
clusively along the [001] direction. Now, however, the 
gold atom hops between pairs of nearest-neighbor O and 
Ti atoms. After ~ 0.8 ps, a configuration is reached 
in which the gold atom is equidistant between the 01 
and 02 (and Til and Ti2) sites before it jumps into an- 
other bridge position between the 02 (green line) and Ti2 
sites. By performing static calculations, we show that Au 
adatoms adsorbed on bridge sites do not induce reduc- 
tion of the substrate (see Fig. [5](B)), thus all the Ti ions 
of the substrate preserve their 4+ oxidation state and a 
very weak oxidation of the Au adatom is observed. This 
is fully confirmed by the dynamical simulations: along 
the trajectory, the Ti-3d orbitals are found to be empty, 
which implies that no localization of charge on substrate 
Ti sites is observed. 

Even at an elevated temperature of 900 K, we do not 
observe diffusion of the Au adatom in the [110] direction. 
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namely from the top of Ob atoms to bridge sites between 
Tisc and Ob atoms, on the timescale of picoseconds. This 
dynamics is consistent with previous findings^^^^ based on 
static or NEB^ calculations which predict a relatively 
high energy barrier, ~ 0.35 eV, for this process to happen 
compared to others. In conclusion, dynamical PBE+U 
simulations demonstrate that Au adatoms diffuse highly 
directionally on the stoichiometric rutile (110) surface. 
They can easily migrate either along the top of the bridg- 
ing oxygen rows of the clean TiO2(110) surface or around 
the area between these rows from one bridging position 
to the next one along the [001] direction. We did not ob- 
serve, on the picosecond timescale, translational motion 
perpendicular to this direction, e.g. from one Ob row to 
a neighboring row via suitable bridging positions. 



Au adatom adsorption on the reduced 
TiO2(110) surface 



D. Au substitutional defects on the TiO2(110) 
surface: AuxTi(i_x)02-(5 



Another reaction channel of gold interacting with ti- 
tania surfaces is via the chemical exchange of Ti atoms. 
We studied the scenario where an Au atom substitutes a 
surface Tisc site, which we call "Au^Vtisc"- The pres- 
ence of such an Au substitutional point defect induces 
a rearrangement of the neighboring atoms, leading to 
the formation of a distorted squared planar "AUO4" unit 
(see Fig. [9]). In this configuration, the Au atom relaxes 
outward by 0.62 A and is found to be coordinated by 
four surface O atoms (at about 2.0 A). The incorpora- 
tion of an Au atom into the titania surface does not yield 
a change in the occupation of the Ti-3(i states: all Ti ions 
preserve their formal oxidation state Ti^+. However, the 
adsorption of an Au adatom into a Vtisc vacancy site 
is strongly exothermic, releasing —6.38 eV (or —6.21 eV 
when using PBE), shown in Table HVl 



Starting with the reduced TiO2(110) surface in the 
presence of an Ob surface oxygen vacancy the Au metal 
adatom is found to adsorb preferentially at the Vo va- 
cancy. Such surface oxygen vacancies result in stable an- 
choring sites for Au adatoms, which bind at about 0.90 A 
above the O vacancy site with two Ti nearest neighbors at 
2.68 A obtained from both PBE+U and PBE. The corre- 
sponding PBE+U (PBE) calculated adsorption energies 
of —1.54 eV (—1.57 eV) are much larger than the bind- 
ing to the stoichiometric surface (see Table HV|) . which 
is consistent (—1.6 to —1.8 eV) with previous GGA cal- 
culations^^. The strong adsorption of an Au atom at 
the Vo site entails a strong charge rearrangement at the 
Au/oxide contact. In the presence of an isolated Vo va- 
cancy, the charge neutrality of the system is maintained 
by the presence of two reduced Ti^+ ions. The bonding 
charge distribution (see Fig. [5](D)) shows that, upon Au 
adsorption at the Vo vacancy site, the charge transfer 
occurs now from the reduced substrate to the supported 
metal atom, thus leaving a reduced surface with a sin- 
gle Ti^+ ion. This indicates that the charge transferred 
from the reduced substrate to the adsorbate comes from 
one of the two Ti^+ ions. As a result this process leads 
to the formation of a negatively charged Au"^" adspecies. 
The analysis of the DOS (see Fig. [5] (D)) now shows 
that the charge transferred from the substrate to the Au 
atom moves toward the latter's half-filled 6s band, which 
turns out to be almost completely filled. In conclusion, 
these calculations not only suggest a greatly increased 
stability of Au adatoms adsorbed onto O vacancies when 
compared to the stoichiometric surface as a reference, but 
also a very different chemical reactivity with respect to 
admolecules in view of their different charge state. The 
resulting ramifications for CO activation will be discussed 
in Sec. IVBl 



This particular defective surface is found to be ex- 
tremely reactive. We have computed the energy required 
to remove one of the oxygen atoms On (n = 1, 2, 3) 
in the surface layer (see Fig. [9] for site labeling). The 
atoms 01 and 03 are Ob atoms while the atom 02 is 
an in-plane oxygen. Our PBE+U values for the 01, 
02, and 03 vacancy formation energies are 1.52, 2.20, 
and 1.36 eV respectively, which is in accord with previ- 
ous PBE calculation^ii. For the stoichiometric, undoped 
TiO2(110) surface, the corresponding PBE+U value of 
Ey is 2.97 eV for an oxygen vacancy, Vo, in the bridg- 
ing row, compared to 01 and 03 here. These values sug- 
gest that substituting a surface Tisc atom with an Au 
atom greatly weakens the binding of surface O atoms, as 
observed in the case of Ce02 surfacesiSSiiii. 



In Fig. E (E), the total DOS and PDOS of the doped 
Au@VTi5c surface is depicted in the presence of an oxy- 
gen vacancy. Here, the missing oxygen is the bridg- 
ing 03 atom, O3 atom (see Fig. [9]), the vacancy for- 
mation energy of which is found to be lowest (1.36 eV) 
when Au substitutes a Ti atom at a Tisc site, which we 
call "Au@VTi5c, V03". Concerning the electronic struc- 
ture, the main difference between the two scenarios, i.e. 
one obtained by substituting a Tisc with an Au atom 
Au@VTi5c and the other generated by substituting a Tisc 
with an Au atom and by removing a surface Ob atom 
Au@VTi5c7 Vo3, is related to the reduction of the ox- 
ide substrate. As we can extract from Fig. [5] (E), in 
presence of a vacancy located in a neighboring bridging 
oxygen row, 03, the excess electron resulting from this 
vacancy transfers to the substrate and a filled state ap- 
pears in the band gap. This gap state stems from a Ti-3d 
orbital, thus reducing one second-layer Ti ion to Ti^+. 
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FIG. 6: Side view (left panel) and top view (central panel) of the diffusion path of the Au adatom on the stoichiometric 
TiO2(110) surface at 900 K. The Au atom, initially adsorbed on top of the bridging oxygen atom 01, diffuses along the Ob 
bridging row by hopping from one oxygen to the next, 01 ^ 02 -^ 03, as visualized by the worm-like trajectory. The 
right panel shows the time evolution of the Au-On (n=l,2,3) bond lengths; the dashed line at 2.00 A refers to the optimized 
equilibrium Au-0 bond length. 




■D 




FIG. 7: Side view (left panel) and top view (central panel) of the diffusion path of the Au adatom on the stoichiometric 
TiO2(110) surface at 900 K. The Au atom, initially adsorbed on a bridge site between the 01 and Til atoms, diffuses by 
hopping in between pairs of nearest-neighbor Ob bridge and Tisc atoms as visualized by the worm-like trajectory. The right 
panel shows the time evolution of the Au-On and Au-Tin (n=l,2) bond lengths; the dashed lines at 2.30 and 2.80 A refer to 
the optimized equilibrium Au-0 and Au-Ti bond lengths, respectively. 



E. Ab initio thermodynamics of defective 
Au/TiO2(110) surfaces 



The effects of temperature and pressure on the rela- 
tive stability of Au/Ti02 metal/oxide surfaces have been 
taken into account by employing the formalism of approx- 
imate ab initio thermodynamics, as sketched in Sec. [Ill 
To this end, we compute the free energies of Au adsorp- 
tion, AGg,ds{T^p) as given by Eq. ([6j), and report them in 
Fig. [To] as a function of the O chemical potential includ- 
ing a conversion to oxygen partial pressures at several 
relevant temperatures. These free energies are measured 
relative to the stoichiometric surface and therefore in- 
clude the free energy cost of creating whatever vacancy 



the gold atom may be associated with. 

As highlighted by the colors in Fig. [TOl it is possible to 
identify four thermodynamically stable phases. The first 
phase, which holds for values of /io > —1.35 eV, cor- 
responds to the scenario where a surface Tisc atom has 
been replaced by an Au adatom, denoted "Au^Vtisc"- 
In oxidative environments, and at the reference con- 
ditions that are traditionally used in most computa- 
tional studies, this structure becomes the thermodynam- 
ically most stable one. The second most stable struc- 
ture thermodynamically is the one obtained by remov- 
ing a surface O atom from a bridging row based on 
the Au@VTi5c structure described above, which is called 
"Au@VTi5c, V03" since the missing Ob atom is the 03 
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t/ps 



FIG. 8: Side view (left panel) and top view (central panel) of the Au adatom adsorbed on top of the surface bridging oxygen 
atom 01, which served as the initial configuration for the migration path depicted in Fig. [6l Right panel: Corresponding time 
evolution of the fractional occupation number of particular Ti-3d orbit als at specific sites, as indicated; populations of about 
and 1 correspond to Ti^^ and Ti^^, respectively. At t = ps, the Au adatom is bound to 01 and the excess charge is fully 
localized in the second subsurface layer at the Ti2 site. 





Tisc atom with an Au adatom, while under Ti-rich con- 
ditions, the Au adatoms are preferentially adsorbed at O 



vacancies. 



FIG. 9: Ball and stick model of an AuxTi(i_x)02 system 
obtained by substituting a surface Tisc atom by an Au adatom 
in side (left) and top (right) views. 



atom (refer to Fig. [9]) We have shown that the PBE+U 
energy required for removing the 03 oxygen atom, i.e. 
the vacancy formation energy Ey^^ from the surface of 
the Au@VTi5c structure amounts to 1.36 eV. This is sig- 
nificantly lower than the required energy of 2.97 eV to 
remove an Ob atom from the ideal, stoichiometric Ti02 
surface. Therefore the TiO2(110) oxide surface becomes 
a better oxidant when doped with gold. 

As demonstrated in Fig. [TOl the Au@VTi5c,Vo3 sur- 
face structure becomes thermodynamically stable for 
-1.79 eV < /io < -1.35 eV. It turns out that 
these two defective surface structures, Au^VxiSc and 
Au@VTi5c, Vo3, are thermodynamically stable in a wide 
range of temperatures T and pressures p that are rel- 
evant for applications in the realm of catalysis. In con- 
trast, the adsorption of Au adatoms on the stoichiometric 
TiO2(110) surface is thermodynamically stable only in a 
quite narrow range of values of the O chemical potential 
of —2.05 eV < /io < —1.79 eV. Finally, the adsorption 
of Au adatoms on O vacancies becomes thermodynami- 
cally stable for values of /io < —2.05 eV. We therefore 
conclude that under 0-rich conditions, the thermody- 
namically most stable structure is the defective surface 
structure Au^VxiSc obtained by substituting a surface 
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FIG. 10: Free energy AGads{T,p) for Au adsorption on stoi- 
chiometric and defective TiO2(110) surfaces as a function of 
the oxygen chemical potential A/io; see text for the four dif- 
ferent phases. Conversion to oxygen partial pressures p (up- 
per axes) has been carried out at constant temperatures of 
T = 300, 600, and 900 K (see text). 
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V. GOLD-PROMOTED TITANIA: 

INTERACTIONS WITH CO PROBE 

MOLECULES 

A. CO molecular adsorption on Au/TiO2(110) 
surfaces 

We have shown that our calculations predict the two 
most stable adsorption sites of an Au adatom adsorbed 
on the stoichiometric TiO2(110) surface, differing ener- 
getically by only ~0.1 eV. The lowest energy configura- 
tion is that with the Au adatom adsorbed onto a bridge 
site between a surface Ob atom and a first-layer Tisc 
atom. A detailed analysis of the electronic structure of 
this configuration shows that all Ti ions belonging to the 
substrate preserve their formal oxidation state Ti^+ . The 
other stable structure consists of an Au adatom sitting 
on top of an Ob atom of the (110) surface. The metal 
adsorption on an Ob top site induces a strong charge re- 
arrangement at the metal/oxide contact and entails the 
reduction of a second-layer Ti ion which becomes for- 
mally Ti3+. 

Using these preferred structures, we now study their 
interaction with a CO admolecule. A CO molecule is 
placed end-on at 2.5 A above the Au adatoms. The ad- 
sorption of the CO molecule on an Au adatom adsorbed 
at an Ob site is found to be strongly exothermic, releas- 
ing -2.66 eV (-2.33 eV using PBE; see Table |IVl). This 
is in accord with earlier calculations^. In this configu- 
ration, the CO molecule is aligned with the Au adatom 
and oriented normal to the surface (see Fig. [11] (F)). The 
computed value of the C-0 bond length is 1.15 A com- 
pared to a value of 1.14 A of an isolated CO molecule 
using the same approach. The PBE+U (PBE) value of 
the distance between the Au atom and the substrate Ob 
atom to which it is bonded is 1.98 A (1.97 A), while the 
Au to C distance is 1.88 A (1.88 A). Analysis of the PDOS 
(see Fig. [11] (F)) reveals that at this Au/Ti02 contact, 
only one Ti^+ ion is present before and after adsorption, 
and that the excess charge populating the substrate oc- 
cupies the same second-layer Ti-3d orbital. Therefore, 
the charge redistribution due to adsorbing a CO molecule 
does not further reduce the oxide support. However, the 
bonding charge analysis (also displayed in Fig. [11] (F)) 
shows that the Au adatom is involved in charge depletion 
(blue areas), while both the C and O atoms accumulate 
the resulting excess charge (red areas). 

We now consider the interaction between a CO 
molecule and the Au atom adsorbed onto the bridge site. 
Again, the mechanism is exothermic by —2.27 eV (or 
-2.17 eV when using PBE instead of PBE+U). In this 
configuration, the CO molecule is bonded to the Au atom 
and the C-0 bond length is 1.15 A as in the previous 
case. Here, however, the CO-Au structure is found to be 
tilted with respect to the surface at an angle of about 60° 
(see Fig.[n](G)). The PBE+U (PBE) values of the dis- 
tances between the Au adatom and the Tisc atom before 
and after the adsorption of molecular CO are 2.79 (2.88) 



and 3.90 A (3.88 A), respectively, while the distances 
between the metal and the O bridging atom before and 
after CO adsorption are 2.30 (2.39) and 1.99 A (2.01 A), 
respectively. The analysis of the PDOS as depicted in 
Fig. [11] (G) shows, in this case, that at the Au/Ti02 
contact, the charge redistribution reduces the oxide sup- 
port, contrary to the previous scenario. We recall that 
when Au is adsorbed onto a bridge site, all the substrate 
Ti ions preserve their formal Ti^+ oxidation state. Upon 
CO adsorption onto the Au atom occupying a bridge po- 
sition, one Ti ion is reduced. The additionally reduced 
Ti^+ ion, which belongs to the first Ti02 layer, is pre- 
cisely the ion that forms a bond with the Au in the bridge 
configuration before CO adsorption. 

The configuration in Fig. [11] (G), where the CO-Au 
structure is tilted by ~ 60° with respect to the substrate, 
is used as an initial condition for a short AIMD run at 
room temperature in order to probe dynamical instabil- 
ities. In less than 0.4 ps, the CO-Au complex bends 
and sits perpendicular to the surface, therefore recover- 
ing the structure obtained upon molecular CO adsorption 
on gold sitting on top of an Ob atom. The difference in 
the binding energies for CO molecule adsorption on the 
top and bridge configurations described above can be at- 
tributed in part to the tilting of the Au-CO complex with 
respect to the surface plane. More important, however, 
is the position of the reduced Ti^+ site. Previous stud- 
^ggi3,79 Q^ charge localization induced by Ob vacancies 
on Ti02 surfaces have shown that excess electrons, which 
are trapped at specific Ti sites, migrate easily by phonon- 
assisted hopping to other Ti atoms, thus exploring dif- 
ferent electronic structure topologies^^. In particular, it 
was found that the most stable sites for charge local- 
ization belong to the second subsurface layer under Tisc 
rows. The topologies where the excess charge is shared 
between surface Tisc atoms and second or third subsur- 
face layer sites below Tisc rows were found to be about 
0.2-0.3 and 0.3-0.4 eV higher in energy, respectively. It 
is therefore conceivable that we would observe a similar 
dynamics for excess charge de- and re-localization. 



B. CO adsorption on Aui/TiO(2-x)(110) and 
Au2/TiO(2-x)(110) surfaces 

The interaction between a CO molecule and a single 
Au adatom adsorbed at a bridging site or on top of an 
Ob atom on the TiO2(110) surface has been investigated 
above. We now concentrate on the interaction between 
CO molecule and Au adatom adsorbed at a bridging 
O vacancy. It has been shown that the metal adatom 
strongly interacts with this reduced oxide support. The 
resulting Au"^" adatom is highly stable, considering its 
binding energy of —1.54 eV. Our calculations predict a 
fairly weak interaction between molecular CO and the 
Au"^" adatom. The CO molecule initially located end- 
on at 2.5 A above the Au'^" species is found to bind with 
about —0.41 eV. This binding energy is close to an order 
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FIG. 11: Electronic structure analyses (based on the PBE+U approach) of (F) a CO molecule adsorbed on a supported Au 
adatom (originally the Au atom was adsorbed at a bridging O atom) on the TiO2(110) surface, (G) a CO molecule adsorbed on 
a supported Au adatom (originally the Au atom was adsorbed between a bridging O atom and a Tisc atom) on the TiO2(110) 
surface, (H) of the Au2 dimer adsorbed onto surface and (I) a CO molecule bonded to an Au2 dimer adsorbed onto a surface 
O vacancy on the TiO2(110) surface. The left panel represents the bonding charge Sp{z) integrated over planes perpendicular 
to the surface and plotted as a function of the height from the surface. The central panel displays the bonding charge Ap{r) 
at an isovalue of ±0.06 |e|/A^ where electron accumulation and depletion are represented by red and blue areas, respectively. 
The right panel shows the total DOS and atom resolved projected DOS (PDOS), as indicated, where energies are with respect 
to the Fermi level, which is marked by a solid vertical line. 
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of magnitude lower than the binding energy of the posi- 
tively charged Au adatom adsorbed on the stoichiometric 
TiO2(110) surface, which binds with -2.66 eV. 

Our calculations show that while positively charged 
Au'^^ ions supported on the TiO2(110) surface are shown 
to activate molecular CO, negatively charged Au"^" ad- 
species, which are incorporated into surface O vacancies, 
interact only weakly with CO adsorbates. This is in 
line with the scenario of Au adsorbed on the Ce02 sub- 
stratei^ii^, where it has been shown that the higher sta- 
bility of Au adatoms adsorbed onto surface is the basis of 
the deactivation mechanism during CO oxidation. The 
oxidation mechanism proposediSS involves three steps: 
the spillover of the CO molecule, the actual oxidation 
via a lattice oxygen atom leading to CO2 desorption, and 
the diffusion of the Au adatom into the newly formed O 
vacancy, which results in negatively charged Au"^" ad- 
species that prevents the adsorption of molecular CO. 
Since the present AIMD simulations have shown that 
Au adatoms diffuse rather easily on the stoichiometric 
TiO2(110) surface, we would expect a deactivation pro- 
cess with titania similar to what was observed on the 
ceria surface. 

Despite the fact that CO only weakly interacts with 
negatively charged Au"^" adspecies, several theoretical 
and experimental studies suggest that Vo vacancies in 
bridging oxygen rows, Ob, are active nucleation sites 
for Aun clusters on the TiO2(110) surfaceiiii^iiii. It is 
therefore expected that nucleation and growth of Au^ 
clusters on the rutile surface is intimately related to the 
presence of surface oxygen vacancies. To get a glimpse, 
we next investigated the limiting case of a gold dimer, 
Au2, adsorbed at an Ob vacancy. The optimized struc- 
ture of an Au2 dimer adsorbed onto a surface bridging 
O vacancy is depicted in Fig. [12] (a). The computed 
value of the adsorption energy of Au2 on the reduced 
TiO2(110) surface, calculated with respect to the iso- 
lated Au2 molecule, is -1.17 eV (-1.19 eV with PBE). 
The adsorption mechanism induces reduction of an addi- 
tional Ti ion, leading to a gold dimer into an O vacancy 
and two second-layer Ti^+ ions (see Fig. [11] (H)). The 
two gold atoms are located 1.20 A above the O vacancy 
site, whereas the distances between the Au atoms and 
the nearest-neighbor Ti atoms are 2.80 A A and the Au- 
Au distance amounts to 2.51 A compared to 2.53 A of 
the isolated gold dimer using the same method. 

Next, a CO molecule is positioned above the Au2 dimer 
adsorbed onto the surface O vacancy. After relaxation, 
the CO molecule reaches the most stable configuration of 
Fig. [12] (b). The computed binding energy is —1.22 eV 
(—1.00 eV using PBE), in good agreement with previ- 
ous studies^^. This value is about half of that obtained 
for adsorption of CO on a single Au adatom adsorbed 
on the stoichiometric surface (—2.66 eV) but distinctly 
lower than the energy calculated on the Au/TiO(2-x) 
surface (—0.40 eV). We therefore find that a cluster as 
small as Au2 nucleated at an O vacancy on the reduced 
TiO2(110) surface favors the formation of stable CO ad- 



sorbates. Specifically, the CO molecule binds to the Au 
atom that is farthest from the O vacancy and the Au-Au 
distance elongates from 2.51 to 2.78 A. The interaction 
between CO and the metal oxide does not induce further 
reduction of the substrate; before and after adsorption, 
two reduced Ti^+ ions are present (see Fig. [11] (I)). In 
Fig- [12 we plot the spin density of the Au2 dimer ad- 
sorbed on the reduced Ti02 surface (a) and the spin den- 
sity after adsorbing a CO molecule to this Au2/TiO(2-x) 
metal/oxide substrate (b). In each panel, we can see that 
both Ti^+ ions belong to subsurface layer sites below Tisc 
rows. 

It is clear from these calculations that a gold "cluster" 
as small as the Au2 dimer adsorbed onto a bridging O 
vacancy leads to a promoted support that strongly binds 
to molecular CO. Furthermore, the step from one to two 
gold atoms greatly changes the electronic properties, thus 
lending support to the general view that the reactivity of 
gold nanoparticles nucleated at O vacancies is strongly 
size dependent. 





FIG. 12: Spin density at ±0.005 |e|/A^ of the Au2 dimer 
adsorbed onto a surface O vacancy on the TiO2(110) surface 
(a), and the same system after coadsorbing a CO molecule 
(b); see text. 



VI. CONCLUSIONS AND OUTLOOK 

We have performed periodic density functional based 
calculations that account for the on-site Coulomb inter- 
action via a Hubbard correction ("GGA+U") on stoi- 
chiometric, reduced, and gold-promoted rutile TiO2(110) 
surfaces. Structure optimizations, ab initio thermody- 
namics calculations, and ab initio molecular dynamics 
simulations have been carried out both using PBE+U 
and plain PBE in order to provide a broad picture on 
the interactions of these surfaces with several catalyti- 
cally important molecules. 

In agreement with plain PBE calculations and experi- 
mental results, we find that H atoms preferentially adsorb 
on surface Ob atoms of the stoichiometric surface. Ad- 
sorption of hydrogen results in its reduction and transfer 
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of close to one electron per H atom to the substrate. Us- 
ing PBE+U, however, the electron transferred by a sin- 
gle H atom into the substrate localizes preferentially at 
second-layer Ti sites (whereas plain PBE yields a delo- 
calized state). Both PBE+U and PBE calculations show 
a decrease in the adsorption energy as a function of cover- 
age and a maximum coverage of about 60-70%. Oxygen 
vacancies are the preferred adsorption sites for H2O dis- 
sociation and similar values for the adsorption energies 
of water are obtained with and without the Hubbard cor- 
rection. Our PBE+U calculations predict the presence 
of two reduced Ti^+ ions before and after the adsorp- 
tion of water, indicating that the adsorption of water 
does not further reduce the support (whereas plain PBE 
again yields delocalized excess electrons). 

Akin to plain PBE calculations, PBE+U predicts a 
weak interaction between CO molecules and the reduced 
Ti02 substrate. The data indicate that CO adsorbs at 
both oxygen vacancies in the bridging surface rows, Ob, 
and fivefold coordinated titania sites in the first layer, 
Tisc. While PBE calculations give similar energy values 
for CO at oxygen vacancies and at Tisc, PBE+U calcula- 
tions show that the adsorption of CO at those Ti^c sites 
which are nearest neighbors of the oxygen vacancies is 
energetically disfavored. Upon CO adsorption on the re- 
duced oxide, the resulting charge redistribution does not 
reduce further the titania substrate. 

Addressing next the interaction of titania with gold, 
both PBE+U and PBE calculations predict two most 
stable adsorption sites for Au adatom adsorption on the 
stoichiometric Ti02 surface. On the one hand, once an 
Au adatom is adsorbed at a bridge site between a surface 
Ob atom and a first-layer Tisc atom, PBE+U suggests 
that all the Ti ions belonging to the substrate preserve 
their formal oxidation state Ti^+ , indicating a very weak 
oxidation of Au. On the other hand, if the Au adatom 
is adsorbed on top of an Ob atom, a net charge transfer 
from the metal to surface, leading to a positively charged 
Au'^^, is observed. In this case, the metal adsorption en- 
tails the reduction of a second-layer Ti ion, which for- 
mally becomes Ti^+. Ab initio molecular dynamics re- 
veals that gold adatoms on stoichiometric TiO2(110) are 
very mobile. They are found to migrate easily along the 
[001] direction, either along the top of bridging oxygen 
rows or around the area between such rows. In the for- 
mer case, we observe an interesting subsurface charge 
delocalization and relocalization dynamics of the excess 



charge. 

Promotion of the (110) rutile surface via adsorption of 
gold atoms or substitution of Ti by Au results in a sys- 
tem with greatly changed electronic properties and thus 
modified reactivities in the realm of heterogeneous catal- 
ysis. This is traced back to the oxidation state of the 
gold adatoms, which is determined by the site at which 
the Au atom is adsorbed as well as by the stoichiometry 
of the substrate. Isolated Au atoms supported by the 
stoichiometric Ti02 surfaces are shown to induce a sig- 
nificant charge redistribution at the metal/oxide contact. 
The calculations, both PBE+U and PBE, show posi- 
tively charged Au'^^ adspecies supported on the stoichio- 
metric surface which activate CO admolecules. In stark 
contrast, the negatively charged Au"^" adspecies, which 
are incorporated into surface Ob vacancies, interact only 
weakly with the same CO molecules. We have shown that 
structures obtained by substituting a first-layer Tisc ion 
with an Au atom weakens the bond of surface O atoms, 
akin to recent observations with ceria. Ab initio thermo- 
dynamics mappings of the surface phase diagram predict 
that the "Au^VxiSc" and "Au^VxiSc, V03" of temper- 
atures and pressures relevant for catalytic applications. 
Finally, we have shown that although a single Au adatom 
bound to a surface O vacancy weakly binds CO, a second 
gold atom adsorbed simultaneously into the O vacancy 
leads to the formation of a gold dimer, Au2, featuring 
distinctly different electronic and binding properties. 

We expect that these insights will be of great help in 
understanding the catalytic activity of gold-promoted ti- 
tania interfaces with liquid water as used in selective ox- 
idation reactions of more complex molecules. 
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